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Abstract. Landau level quantization in graphene reflects the Dirac nature of its 
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holes and this leads to characteristic behaviour of the magneto-optical conductivity as 
a function of frequency f2 for various values of the chemical potential /i. Particular 
attention is paid to the optical spectral weight under various absorption peaks and 
its redistribution as ^ is varied. We also provide results for magnetic field B as well 
as chemical potential sweeps at selected fixed frequencies which can be particularly 
useful for possible measurements in graphene. Both diagonal and Hall conductivity 
are considered. 
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1. Introduction 

Recent experimental studies pj, i2j of the dynamics of electrons and holes in graphene 
(a single atomic layer of graphite [3]) have revealed unusual behaviour related to the 
Dirac nature of its quasiparticles. Two dimensional graphene has a honeycomb lattice 
structure with two atoms per unit cell. Its band structure consists of two inequivalent 
pairs of cones with apex at the Brillouin zone corners. For zero chemical potential the 
lower energy cones are completely filled and the upper empty. In a graphene device an 
applied gate voltage can be used to introduce electrons in the upper band or, by voltage 
reversal, holes in the lower band (cones). 

The quasiparticles in graphene obey Dirac [H [5] rather than the Schrodinger 
equation and this has profound implications for their dynamics. The unconvetional 
quantum Hall effect was expected theoretically P, [71 El El and recently observed [H [2] 
to have half (divided by spin and valley degeneracy) rather than integer filling factors. 
The predicted phase shift of vr in the de Haas van Alphen [HI [121 IE] and Shubnikov de 
Haas [Hj oscillations was also seen experimentally [H [2] . Another feature related to the 
Dirac-like character of the carriers in graphene seen in the dc measurements is a finite 
effective cyclotron mass for the massless Dirac quasiparticles which varies as the square 
root of the number of carriers pTl [T^ [U [2]. 

In this paper we consider the magneto-optical conductivity of graphene. Work 
without a magnetic field includes the calculations of Ando, Zheng and Suzuura [15] who 
considered the effect on the frequency dependent conductivity of short- and long-range 
scatterers in a self consistent Born approximation. More recent work [16] describes 
several anomalous properties of the microwave conductivity of graphene with, as well 
as without magnetic field. These properties are directly related to the Dirac nature 
of the quasiparticles. Several analytic formulas for the longitudinal as well as Hall 
ac conductivity are given in the paper [8]. They also present extensive results for 
dc properties and preliminary data on the real part of a^xi^) vs Q in the optical 
region. Another extensive work by Peres et al. p] on the ac conductivity in graphene 
treats localized impurities in a self-consistent fashion as well as extended edge and grain 
boundaries including also effects of electron-electron interactions and self-doping. 

In this paper we follow most closely Ref. [8] which we extend in several directions. 
An aim is to provide simpler analytic formulae which should prove useful in the analysis 
of experiment and check on their accuracy. Another is to consider magnetic field as 
well as chemical potential sweeps possible in graphene field effect transistor devices 
[H [21 [31 [17]. Recent work by Li et al. in organic metals [T8| has demonstrated 
that the ac measurements are also possible in such devices. Consistent with our 
aim, we consider impurities only in simplified scattering rate model and neglect real 
part renormalizations although these could easily be included if wished. While the 
renormalization effects beyond those included in this simplified model might become 
important for the interpretation of future experiment, we point out that so far the free 
quasiparticle model with associated transport lifetime has been remarkably successful 
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in understanding the dc results of Refs. [H [2]. 

In Sec. [2] we relate the magneto-optical conductivity tensor to the Dirac fermionic 
Green's function through Kubo formula. The general formulas obtained can be greatly 
simplified and closed form expressions are obtained in two cases. In general the fermionic 
self-energy can depend on energy and Landau level index as well as temperature and 
value of the external magnetic field B. Under the assumption that variations with 
Landau level index n can be neglected, the sum over transitions between neighboring 
Landau levels can be carried out explicitly and a closed form expression is obtained for 
conductivity in terms of digamma function. A single integral over an internal frequency 
remains. The expression obtained is suitable for calculations of the ac conductivity 
for any value of temperature, chemical potential and magnetic field. Its microwave 
frequency limit has been used to describe properties of graphene in Ref. [Hj. In the 
case when energy dependence of the fermionic self-energy can be neglected, the internal 
integration over energy can be done and what remains is the sum over the Landau 
level index n of Lorentzian forms multiplied by thermal factors and algebraic weighting 
factors. The weighting factors depend on the Landau level energies as well as the 
excitonic gap (see e.g. Refs. [191 EQl El]) should one wish to include this possibility. 
When we compare numerical results obtained from Lorentzian model and from the 
previous more general expressions in the limit of constant scattering rate, we find good 
quantitative agreement between the two. This provides support for the analysis of 
experimental data |22]. In Sec. 12.31 we consider the low field limit for the Lorentzian 
model derived in Sec. 12. 21 and establish its correspondence with previously known results. 
In Sec. I3.1l we present the numerical results for the real part of the diagonal conductivity 
as a function of photon energy Q for fixed value of magnetic field and various values of 
chemical potential. We also provide results for fixed photon energy while sweeping either 
chemical potential or external magnetic field which should prove useful in comparing 
with experiments. The effect of opening of an excitonic gap on the absorption lines 
is described. In Sec. 13.21 a discussion of optical spectral weight redistribution by the 
magnetic field is given. Section 13.31 is structured in parallel to Sec. 13.11 but deals with 
the absorptive part of the off-diagonal Hall magneto-optical conductivity. Discussion 
and conclusions are given in Sec. HI Some of the algebra needed in this work is found 
in an 



Appendix A 



2. Analytic expressions for optical conductivity 

The optical conductivity tensor is calculated using the Kubo formula 

^ii(^) = — ' (1) 

where 11^ is the retarded current- current correlation function which in the bubble 
approximation is given by 

n,,(l]+^0) = 6^4 r dudu'^^^^^^^ f ^tr [YAiu;, kh^A{u;\ k)] , z = 1, 2(2) 
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where np{u) is the Fermi distribution function uf^u) = l/[exp((c<j — n)/T) + 1], tr not 
only takes care of the 4x47 matrices [7*^ = 0-3^ (as, ia2, —iai)], but includes also the 
summation over flavour (spin) index. Here 



n=0 



(7°M„ + A)/i(k) + /2(k) ^ (7°M„ - A)/i(k) - /2(k) 



is the spectral function (decomposed over Landau levels) associated with the translation 
invariant part of the Dirac fermion Green's function in an external magnetic field B 
applied perpendicular to the plane along the positive z direction (see e.g. Refs. [2T|[TT]). 
In Eq. ([3D 



/i(k) 



P-Ln ( -j — 7-7 ] — P+Ln-l 



eB - — V e5 



.2\ 1 



2ck2 

JeBl 



(4) 



with P± = (1 ± z7^7^sgn(ei?))/2 being projectors and L^{z) the generalized Laguerre 
polynomials. By definition, Ln{z) = L^{z) and L°]^(z) = 0. The energies of the 
relativistic Landau levels in Eq. ([3]) are 



En = ±Mn, M„ = ^^A^ + 2nv^^\eB\/c, (5) 

where the energy scale associated with the magnetic field expressed in the units of 
temperature reads 

^ ^ ^^-^ = 8.85 X 10~'K\Um/s)BiT), (6) 

where vp is the Fermi velocity in graphene given in m/s and the field B is given in 
Tesla. In the following we set ^ = fc^ = 1, and in some places e = c = 1, unless stated 
explicitly otherwise. For the numerical calculations we use the value vp ~ lO^m/s [HIS] 
which leads to the relationship eB (8.85 x 10^K^)i?(T). We consider relatively low 
fields B < 17T, where spin splitting is unresolved [23], so that we can assume that the 
above mentioned summation over flavor index simply gives Nj = 2 in all expressions 
below. The Landau level energies ([5]) include also an excitonic gap A. The physical 
meaning of this singlet excitonic gap is directly related to the electron density imbalance 
between the A and B sublattices of the bi-particle hexagonal lattice of graphene [T9l [20] 
and there are strong indications [2l| that it was indeed observed in recent experiments 
[23[ [25] . We will see here that optical measurements done on graphene can be very 
useful in investigations of the excitonic gap. 

Finally, the scattering rate r„(u;) is expressed via the retarded fermion self-energy, 
r„(a;) = — ImS^(ci;) which in general depends on the energy, temperature, field and the 
Landau levels index n. This self-energy, which in general has also a real part, has to 
be determined self-consistently from the Schwinger-Dyson equation. This equation can 
be solved analytically [20] and numerically as was done in Refs. [6l [9] . In our paper we 
consider Tn{uj) as a phenomenological parameter for the two cases: (i) T{uj) = Tn{uj) 
is independent of the Landau level index n and (ii) r„ = Tn{uj) is independent of the 
energy uj. Under these assumptions the optical conductivity can be studied analytically. 
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2.1. Frequency dependent scattering rate 

Assuming that r„(ci;) is independent of the Landau level index, i.e. T{uj) = r„(c<j), one 
can calculate the sum over Landau levels and express the diagonal ac conductivity in 
the closed form [8] 

e^N ( IB 

ReaU^) = ^ j duiuFiu) - np{uj')]Re ^^—^ [^i{-B) - ^^(-i?)] (7) 

+ + - E,{-B) - E,{+B)] i, i f'^^j^'^^' ^ + ^ cD', F ^ F') 

Here tp is the digamma function, we denoted B = Vp\eB\/c, and also included the 
renormalization of energy caused by the real part of self-energy uj{uj) = u — ReS((D), 
uj' = u{uj + Q) , r = F(u)), F' = and introduced the following short-hand notations 

(cu' + iF')(cD + iF) - A2 



[±B) = E,{u,u',r,r',±B) 

[±B) = E,{u,Cu',T,T', ±B) 



[u-u' + i{T - F')] p + cD' -I- i{T + F')] ± 2B ' 
{u' -tr'){u + ir)-A^ 



[Co-u' + i{r + r')][u + u' + i{r -r')]±2B' ' ' 

The advantage of Eq. ([7]) is that the ip function contains the contribution to conductivity 
from all transitions between neighboring Landau levels. An infinite number of these 
transitions need to be taken into account when the limit i? ^ is considered, so that 
the zero field limit is easily treatable [8] on the base of Eq. (I?]). Another important 
feature of Eq. (^^) is that we kept the frequency dependent impurity scattering rate F(ci;) 
which allows us to investigate its influence on the shape of the Drude peak [16] . 

Although in our work we will mostly use Eq. (JTj) for the numerical computations of 
the diagonal conductivity, there is a possibility to derive simple approximate expressions 
for the diagonal and Hall conductivities which turn out to be very useful when one is 
interested in the resonance peaks of these conductivities in the infrared region. 

2.2. Landau level index dependent scattering rate and magneto- optical Lorentzian 
model 

For analyzing experimental data it is very useful to have a magneto-optical Lorentzian 
model for the complex conductivity, cr±{Q) = axx{^) ± icTxyi^) due to inter- and 
intraband Landau level transitions [26j. For the case of Dirac fermions it is considered 
in 



Appendix A, whereof we obtain the complex diagonal conductivity 



5^ I ( 1 - M M +J ~ np{Mn+i)] + [n^(-M„+i) - n^(-M„)]) 

1 1 \ 1 



Mn - Mn+l + + liVn + F„+i) M„ - M„+i liVn + F„+i) J M„+i - M, 

/ A2 \ 1 
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and the complex Hall conductivity 

a,y{Q) = / V ([ni.(M„) - ni.(M„+i)] - [ni.(-M„+i) - n^(-M„)]) 
Znc ^ — ' 

n=0 



X 



A 



2 



MnMn+lJ Mn+l-Mn 
1 



+ 



M„ - M„+i + + liVn + r„+i) M„ - M„+i - - z(r„ + r„+i; 

A2 \ 1 

1 + 



Mn + M„+i + + z(r„, + r„,+i) ^ M„ + M„+i - - z(r„ + r„+i) ) } " 

In deriving Eqs. ([9]), ( ITOi) we assumed that the impurity scattering rate is independent 
of the energy cu, but kept its dependence on the Landau level index n. This assumption 
allows us to eliminate the integration over uj which is present in Eq. ((Tj). To preserve 
the Landau index dependence, the sum over transitions between neighboring Landau 
levels is retained in Eqs. ([9]), ( fTOl) . 

Based on Eqs. ([9]), ( ITOi) one can easily write down separate expressions for Re axxiyi)) 
\m.axx{^)-, R'e(Ta;j^(f2), \maxy{^) and verify that diagonal and off-diagonal conductivities 
satisfy Kramers-Kronig relations. A big advantage of Eqs. (j9]), (|TOl) is that they are more 
suitable for numerical computations and it is sufficient to include only a few terms in 
the sum even in relatively low magnetic field. Also they are useful for the description of 
the resonance peaks when the Landau level index dependence is more important than 
the energy dependence of the scattering rate which is included in Eq. ([7]). 

All figures in this paper were computed for F = const, so that for calculation of the 
diagonal conductivity we used Eq. ([7]) and checked that it gives almost identical results 
to Eq. ([9]). Since the corresponding full expression for the optical Hall conductivity 
derived in Ref. [8] is rather difficult for numerical calculations, in the present paper we 
compute all results for the Hall conductivity using Eq. ( fTOl) . In the case of Landau level 
index independent width, r„ = const, Eqs. ([9]), (|TOl) acquire an even simpler form 



\ [npi-M^) - n^(M„+i)] + [n^(-M„+i) - n^(M„)] 1 



and 

axy{n) = ^ V {[npiMn) - UFiMn+i)] - [npi-Mn+i) - npi-Mn)]) (12) 

TTC ^ — ^ 

n=0 
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One can see that the conductivity a^xi^jf^) is even function of fi while axy{^,fi) is an 
odd one. 



2.3. Low field limit of the magneto- optical Lorentzian model 

Now we check that Eqs. fllll) . (fT2|) reproduce correctly the limit B ^ 0. Introducing the 
continuum variable oo instead of M„ given by Eq. (JSj) and replacing the sum over n by 
the integral, we obtain 



2ie^{n + 2iT) 



h 



duj- 



UJ 



duo 



_{Sl + 2iVf uj 
(fi + 2iYf - 4u;2 



dni 



dtu 



(13) 



and 

e^vleB p fdnpicu) , dnp{-co) 

' tu^-A^ 1 . a;2 + A2 



(14) 



X 



+ 



(f] + 2zr)2 4cj2 - (f] + 2ir)2 



where we restored Planck constant h = 2'Kfi in the overall prefactor. Here the terms 
with the factor (cj^ — A^)/^;^ are intraband and the terms containing the factor 
(tu^ + A^)/^^ are interband. The expressions (fT3l) and are obtained under the 
condition h\eB\v'^p / c <^ F. The intraband term of Eq. (|T3l) can be written in the 
familiar Drude form 



Drude 

XX 



2e' 



du 



dnp{uj) 



{uj^ - A^)d{u^ - A' 



(15) 



h J -oo \ duj J 2T — iil \uj\ 
which for A = reduces to Eq. (2) of Ref . [16] . The whole expression (fT3|) which includes 
the interband term, also reduces to the other limiting cases considered in Refs. [TB1[27|. 
In particular, an unusual feature of graphene is that in the high frequency limit the 
interband contribution is a constant, 

a,,(fi)~— , n>/i,A,T. (16) 
2n 

Here to rely on the linearized Dirac approximation we assumed that Q is still well below 
a large band edge. 

The real part of the Hall conductivity (HM for Q = A = T = reduces to the 
expression 

e'^VpeBsgnfi 



axy(^ = 0) 



h\eB\vjp/c < r < 



:i7) 
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which is in agreement with Eq. (4.3) of Ref. [5]. On the other hand, in the high frequency 
hmit Eq. (|T^ gives 

A" 



tanh h tanh ^ 



n 



oo. 



1^ 



2T 2T 

This behavior also follows from Eq. flTU]) which is valid in an arbitrary field B. 
Interestingly, expression (fTSjl is sensitive to the relationship between and A and 
this feature can be used for the gap detection [see the discussion of Figs. [Hand [9] below]. 
Using Eqs. f|T6l) and f|T8|) we obtain the weak field optical Hall resistivity 



IQv^pB 



tanh 



A , u 

+ tanh ■ 



A' 



2T 2T 

Accordingly, the high frequency optical Hall coefficient for A = and T 



sgn(/i) 



(^>/i,A,T. (19) 


(20) 



B necn^ 

contains the information on the value of the Fermi velocity in graphene. 

Finally, the imaginary part of cr^y(f2) which follows from Eq. f[T^ is given by the 
expression 

oo 

ie^leBTQ 



Imcr^y(fi) 



TTC 



dnF{uj) ^ dup 



duo 



du 



A 
1 



+ 



1 



(21) 



(^2 _ 4r2)2 + i6f]2r2 [Ay^ - Q? + AT^ f + IQ9?T^ 
which in high frequency limit becomes 

The last equation shows that in this limit real part of cr^y given by Eq. (ITS!) is the leading 
term. 



3. Results for optical conductivity 
3.1. Diagonal conductivity 

In Fig. [1] we show the results based on a numerical evaluation of the full equation ([7]) for 
the real part of the longitudinal conductivity Kea^xi^) in units e'^/h as a function 
of frequency in units of cm~^. Except for the long dashed (red) curve which was 
obtained in the limit of vanishing external magnetic field (namely B = 10~^T) and 
is included for comparison, the other three curves are for B = IT. They differ in 
value of the chemical potential /i. In all cases the temperature T = lOK, the impurity 
scattering F = 15K, and the excitonic gap A = 0. For reference in scrutinizing the 
curves, the frequency of the n = 1 Landau level, Qi = Mi (A = 0) = 294cm~^[423 K], 
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^2 = M2(A = 0) = 415.8cm-i[598K] and = M^iA = 0) = 509.2cm-i[733K], so 
that for the dash-dotted (black) curve fi falls below the energy of the n = 1 level (see 
leftside of Fig.[2](c))lil, for the solid (blue) curve it falls between n = 1 and n = 2 levels 
(see middle of Fig. [2] (c)), and for the short dashed (green) it is between n = 2 and 
n = 3 (right side of Fig. El (c)). The energies of the peaks are Mi, Mi + M2, M2 + M3, 
etc. and M2 — Mi, M3 — M2, etc. Note that the position and the intensity of the last 
two peaks in Fig. [1] (largest frequency Q) remains the same for all chosen values of the 
chemical potential /x. When /x falls between M2 and M3 [short dashed (green) curve] 



^^=50K B=0.0001T 




200 400 600 800 1000 1200 

n[cm"''] 



Figure 1. (Colour online) Real part of the longitudinal conductivity, Keaxxi^) in 
units of e^/h vs frequency fl in cm~^ for temperature T = lOK, scattering rate 
r = 15K. Long dashed, the chemical potential /i — 50K and the magnetic field 
B = IQ-^T, dash-dotted fi = 50K and B = IT, sohd ^ = 510K and B = IT, short 
dashed = 660K and B = IT. 

the intensity in the third last peak has dropped to half the value it has in solid (blue) 
curve, while the fourth last peak has merged into a low intensity background, as has the 
fifth last, which is seen only in the dash-dotted curve. Also a new peak has appeared at 
M3 — M2 which was not present in the dash-dotted curve. Similarly for fi between Mi 
and M2 [solid (blue) curve] the intensity of the fourth highest energy peak has dropped 
to half the intensity it has in the dash-dotted curve and the peak at Mi is entirely 
missing, having merged into the low intensity background. A new peak has appeared 
at M2 — Ml. Finally when /i is below Mi (dash-dotted curve) there is no peak below 
Ml and the line at Mi always has full intensity. Whatever the value of fi this line will 

X We recall that the conversion rule from the frequency in cm^^ to the energy in Kelvins is 
n[K] = 1.4387K • cmr2[cm-i]. 
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never be seen to half its intensity and this is the hallmark of the Dirac nature of the 
quasiparticles . A schematic which helps us understand the behaviour of the absorption 
hues that we have just described is given in Fig. [2] (c). On the left of the figure we 
show the energies of the Landau levels En = of Eq. ([5]) as solid (red) dots along 

with their quantum numbers n = 0,l,2,.... The Dirac cones which come in pairs, with 
the positive and negative energies and exist at K and K' points in graphene Brillouin 
zone are also shown. The three values of chemical potential considered in Fig. [2] (c) 
are shown in heavy solid (violet) horizontal lines. The possible optical transitions in 
each case are indicated as vertical arrows and connect levels n to n ± 1 only. Moving 
from left to right we see first a single transition from Eq to Ei = Mi, then a pair of 
interband from Ei = —Mi to E2 = M2 and E2 = —M2 to Ei = Mi followed by another 
pair from E2 = -M2 to E3 = M3 and E3 = -M3 to E2 = M2 etc. For the middle 
set of lines the first is an intraband transition from Ei = Mi to E2 = M2 followed by 
a single interband from Ei = —Mi to E2 = M2 and then a pair from E2 = —M2 and 
E3 = M3 and E3 = — M3 to E2 = M2 etc. Finally in the set of transitions at the right 
of the figure there is an intraband from E2 = M2 to E-^ = M3 followed by a single 
interband from E2 = —M2 to E^ = M3 and a pair from E^ = — M3 to E4 = M4 and 
E4 = — M4 to E3 = M3 etc. This is precisely the pattern we have seen in Fig. [H We 
note one more aspect of the anomalous line at the far left of the figure for the transition 
form n = to n = 1 Landau level. It is the only line which cannot be unambiguously 
assigned to inter or intraband because the sate at n = falls at the apex of the Dirac 
cones, where positive and negative energy cones meet and hence they share this state 
equally. This is further illustrated in Figs. [2] (a) and (b), where two sets of cones at K 
and K' are shown separately for the case of A = in frame (a) and finite excitonic gap 
in frame (b). In this second case we see clearly that points K and K' react differently 
under a finite magnetic field. For the cone on the left n = level has moved to energy 
Eq = —A and for the cone on the right it has moved to Eq = A. Note that for the 
value of chemical potential shown as a solid horizontal (violet) line, the Eq = —A to Ei 
transition (vertical arrow) on the left is now definitely interband and the = A to Ei 
on the right is intraband. The ambiguity in designation of the n = level present in 
the top frame (a) of Fig. [2] is lifted when the gap becomes finite. 

The rather complicated pattern of behaviour just described can be understood 
simply from the mathematics of the previous section in the limit F — > and T ^ 0. 
Taking A = 0, yU > 0, and 1] > one obtains from Eq. (1A.12I) 




6{Mn + Mn+i - n) 
Mn + Mn+l 



(23) 
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Figure 2. (Colour online) (a) Schematic representation of the two pairs of Dirac cones 
with apex at points K (left) and K' (right) in graphene Brillouin zone. The energies 
of the Landau levels are shown for index n = 0, ... 4 as solid (red) circles for both 
positive and negative Dirac cones. The transition from n = to n = 1 across the 
chemical potential shown as a thick horizontal (violet) line are for the case A = 0. (b) 
Same as (a) but now there is a finite excitonic gap A. (c) The pair of cones at points 
K and K' in the Brillouin zone (see panels (a) and (b) ) are combined. Vertical arrows 
show allowed optical transitions between Landau levels for three values of the chemical 
potential. At left are the transitions when n is between n = and n = 1, in the middle 
between n = 1 and n — 2, and the right between n = 2 and n = 3. The first line in the 
left is different in that it is both inter (between two separate cones) and intra (within 
a given cone). The first line in all other series in intra and the shortest inter appears 
only once, while all others appear twice. Note that from Eq. (jlip axx{^,fJ') is even in 
^, while from Eq. (fT2|l axy{^, fJ^) is odd in /i. 
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For /i E]Miy, Miy+i[ the T = thermal factor 
2 - npiMn) - npiMn+i) = { 1 



for n < N, 

1 for n = N, 

2 for n > N, 



(24) 



while 



1 for n = N, 
for n ^ A^. 



(25) 



The line in lieaxxi^) with the frequency M„ + M„+i occurs only for n > N with the 
others missing. Further the line with n = N has half of the weight it would have if 
it occurred in another case, namely, 1/(M„ + Mn+i), while lines for n > have full 
weight 2/{Mn + M^+i). Note that the lines with the frequency equal to the difference 
in Landau level energies have weight l/(M„+i — Mn). These lines are always present 
except for the case when /i falls below Mi when the single line at Mi has weight 2/Mi. 
In summary, we have seen in the above discussion that as /i moves through higher and 
higher values of Mjv the lines below disappear into the background with a new line 
appearing at M^v+i — Mj^f. Further the line at n = loses half its intensity while the 
others remain unaltered except for the special case when chemical potential falls below 
the Landau level Mi. 

The pattern would be quite different if instead of M„ ~ ^/n the Landau level 
quantization was Schrodinger-like M„ = uJc{n + 1/2) with being the cyclotron 
frequency. In this case the position of the line corresponding to the difference in Landau 
level energies (intraband) never shifts in energy. Further, the lines corresponding to the 
sum of the Landau level energies (interband) fall at regular energies intervals 2u}c{n + l), 
namely 2uJc, ^oJc-, BcUc, .... Furthermore as increases through the energies of the various 
Landau levels all lines half their intensity before fading into the background and for /i 
below the first level there is no line at uOf.. There is also no ambiguity about whether a 
line is intra or interband as no level is shared between upper and lower cone. 

In field effect devices the chemical potential in a graphene sheet can be changed 
by applying a gate voltage and this may be an ideal way to observe the effects just 
described. However, there is an alternative way to see the same effects. For any fixed 
value of chemical potential, the external magnetic field can be selected in such a way 
that /i falls below the first Landau level or in between n = 1 and n = 2 etc. For 
appropriate choices of B [28] the curves for Re(Ta,a,(n) can be made to behave exactly 
as in Fig. [H To see this it is important that Vl be divided by a/ heBvp / c, so that the 
lines remain fixed in normalized frequency and the vertical scale be divided by the same 
factor and multiplied by F to keep the dimensions the same. When this is done, the 
same pattern as seen in Fig. [1] emerges for this configuration corresponding to fixed /i 
with several well chosen values of B. 

We have found that the curves of Fig. [1] change very little as the chemical potential 
is varied within the limited range M^r to M^r+i. We have also verified that for the 
parameters used here T = lOK, F = 15K the crossover from half intensity in the line 
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Mjsi + Mj\i^i to near zero, i.e. merging into the background occurs rather abruptly in an 
energy range set by T and/or F. This is illustrated in Fig. [3] where we consider how the 

^i=^OOK B=1T 

fi=4^0K B=1T 

,[/=420K B=1T 

,[/=430K 8=1 T 



M 

200 400 600 800 

n[cm"^] 

Figure 3. (Colour online) Real part of the longitudinal conductivity, Re(Jxx{^) in 
units of e^/ft. vs frequency fl in cni^^ for temperature T = lOK, scattering rate T = 15K 
and magnetic field B = IT for 4 values of chemical potential. Long dashed ^ — lOOK, 
dash-dotted ^l = 410K, solid fi = 420K, short dashed fi = 430K. 

complete disappearance of the peak at i7 = fii = 294cm~^ (420K) and the depletion of 
the next higher peak at 710cm^^ towards half its initial value proceeds as the chemical 
potential fi crosses through the energy of the n = 1 Landau level [Qi = Mi (A = 0)]. The 
long dashed curve (red) is for fi = lOOK chosen to be well away from the crossover point 
of 420K and is shown for comparison. The dash-dotted curve (black) is for fi = 410K 
slightly below the crossover point, the solid curve (blue) is for = 420K just at the 
crossover energy and the short dashed curve (green) is for /i = 430K which spans ±10 K 
on either side of Qi which is much less than the level width of 2F = 30 K. As n increases 
through 410K we note the growth of the peak at = 122cm~^, the depletion towards 
zero value of the peak at 294cm^^ and towards half its initial value of the next higher 
peak at 710cm~^. The complete transfer of spectral weight between the various peaks is 
completed for a rather small region of chemical potential variation about Qi = 294cm~^ 
(420K) of order 30 K (not shown in the figure). After the crossover is complete, the 
pattern of the spectral weight distribution will remain unchanged until the chemical 
potential becomes close to the energy of the next Landau level. We would like to stress 
that although Fig. [3] is plotted using the full expression ^ , we have verified that the 
results obtained using Eqs. ( 1A.12I) and ([9]) are practically identical. 

Returning to Eqs. (123|) to (!25|) we note that the optical spectral weight under a given 




Magneto- optical conductivity in Graphene 



14 



inter- or intraband line varies as the square root of B (see Eq. f l32l) and the discussion 
associated with it). This has been verified in the recent experiment of Sadowski et al. 
[29j on ultrathin epitaxial graphite samples [30l [31] . 



H=^5QK A=OK 
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Figure 4. (Colour online) Real part of the longitudinal conductivity, Tleaxxi^) in 
units of e^/h vs frequency O in cm"-'^ for temperature T = lOK, F = 15K, B = IT 
and chemical potential /i = 150K for 4 values of the excitonic gap A. Long dashed 
A = OK, dash-dotted A = lOOK, solid A = 150K, short dashed A = 200K. 

The effect of an excitonic gap A on the optical spectral weight distribution in a one 
Tesla magnetic field is illustrated in Fig. |H Here as in the previous figures F = 15K and 
T = lOK. The chemical potential is set at 150K in all cases. The long dashed (red), 
dash-dotted (black), solid (blue) and short dashed (green) curves are for A = 0, lOOK, 
150K and 200K respectively. The curve for A = is for reference. We note that for finite 
A = lOOK the line at = Mi(A = 0) = 294cm"^ splits into two peaks [dash-dotted 
(black) curve]. The lower peak is at energy Mi (A) —Mq{A) = a/^^^ + A^ — A, while the 
upper peak is at Mi(A) -|-Mo(A) = a/^^? + A^ + A. These two transitions for that value 
of chemical potential are illustrated in Fig. [2] middle frame (b) (see arrows). Additional 
transitions not shown in this frame can, of course, occur but these will have higher 
energy. For A = 150K = fi the n = state shown on the right hand cone of Fig. [2] (b) 
is occupied with a probability 1/2 at T = 0, F = and a transition from n = 1 lower 
cone to n = upper cone is possible as is from n = upper cone to n = 1 of the same 
cone and each must be weighted by a factor 1/2. The first has energy ^nj + A2 + A, 
while the second has energy ^nj + A2 - A. There is an additional transition coming 
from the second cone on the left hand side of the figure of energy v/fif + A2 + A. Thus 
the spectral intensity of the lower energy line in the solid curve (blue) is lower than that 
of the higher energy line by a factor of 3. Finally for A = 200K (A > /i) the energy 
of the lowest transition possible in both left and right side cones are the same equal to 
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^y^lf + A'^ + A and so there is only one line in the dashed (green) curve. This pattern 
of behaviour should allow one to measure the occurrence of an excitonic gap catalyzed 
by the magnetic field. We note that the splitted peak in Fig. H] with finite A acquire a 
slight asymmetry because of the factor 1 ± A^/M„M„_(_i of Eq. (]A.12|) which is different 
for inter and intraband transitions. 

So far we have shown the results for fixed values of chemical potential and magnetic 
field as a function of photon energy. It is also of interest to fix the photon frequency and 
sweep either chemical potential (/i) or field (B) [32j. We begin with the first of these two 
possibihties and this is illustrated in Fig. [51 Four well chosen frequencies are selected, 
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Figure 5. (Colour online) Real part of the longitudinal conductivity, Keaxxi^) in 
units of /h vs the chemical potential /i in K for T = lOK, T = 15K, B = IT and 
A = 0. Four frequencies are considered, long dashed — 122cm~^, dash-dotted 
n = 294cm-\ solid n = 710cm-\ short dashed n = 0.95^1 = 280cm-i. The first 
three frequencies correspond tofi = fi2 — i^i, f^i and CI2 + f^i, respectively, for the 
parameters used. 

viz. Q = Q2 ~ ^1 = 122cm^^ long dashed curve (red), = = 294cm^^ dash-dotted 
curve (black), Q = Q2 + ^1 = 710cm~^ solid curve (blue), and Q = 0.95Qi = 280cm~^ 
short dashed curve (green). Taking these in order, we note that the long dashed curve 
(red) is near zero until fi reaches the value of 420K at which point it increases rapidly 
reaching a plateau at Reaxx{^) ~ 17.5e^//i [the height being set by the height of the 
peak in the solid (blue) curve in Fig. [1] after which it drops rapidly as /i goes through 
595K. It does not drop down all the way to zero however as it continues to sample 
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parts of the peaks at the difference in frequencies M„+i — M„ until these move below 
the sampling optical frequency set at 122cm~^. The next optical frequency chosen at 
f2 = f2i short-long dashed (black) curve begins by sampling the equivalent curve in 
Fig. [1] where the peak has height ~ lAe^/h. However when /i crosses the energy of 
the first Landau level at 420K this peak disappears and the curve drops to zero. No 
other level crosses this frequency again. The short dashed curve is for Q = 0.95f2i. 
It shows very much the same behaviour as the previous case, but its plateau height 
is a little smaller because the peak (dash-dotted curve of Fig. [1]) is sampled not at its 
center, but rather slightly below its maximum value. The final solid (blue) curve is for 
Q = Q2 + ^1- In this case a first plateau is seen at small fi with height ~ b.be^/h 
which is the height of the second peak in the dash-dotted (black) curve of Fig. [H It 
drops to half of this value as VL crosses 420K and then near zero as Vl crosses 595K. 
All features of these curves can be traced to the corresponding behaviour of the curves 
of Fig. [H An additional feature of these curves in now described. Note that as the 
chemical potential gets small the black (dash-dotted) curve increases slightly, while the 
green (dotted) curve drops. As we have described in connection with Fig. [3] the line 
shapes do depend on the value of /i if it falls within order F and/or T of a Landau level 
energy M„. For near zero, the first peak height in the black (dash-dotted) curve of 
Fig. [1] increases by ~ 5% and its width is also slightly narrowed. This leads to a slight 
increase in peak height at = i7i = 294cm~^ monitored in the black (dash-dotted) 
curve of Fig. [5] and a reduction in the green (dotted) curve which monitors the height 
of the curve slightly off the peak at Q = 0.95f2i = 280cm~^. 

Another set of useful curves when considering possible experimental configurations 
is to fix the frequency of the light as well as the chemical potential and sweep the 
magnetic field. Results are shown in Fig. [6] for three cases VL = 200cm~^ long 
dashed (red) curve, dash-dotted (black) curve for Q = 400cm~^ and solid (blue) for 
Q = 800cm~^. In all cases T = lOK, F = 15K, /i = 50K and the excitonic gap A is 
taken to be zero. On the horizontal axis we plotted 1/^B{T), the inverse square root of 
the magnetic field in Tesla, so that fields below one tesla fall above one on this scale. The 
pattern of oscillations is perhaps more complex than in previous curves, but can easily be 
traced out from a knowledge of such curves for Re a^xi^) at different values of B. As an 
example consider the case f2 = 200cm~^. The value of Iieaxx{^)h/e'^ at 1/ a/-B(T) = 1 
is just that of the dash-dotted (black) curve of Fig. [1] read at Q = 200cm~^. As B is 
decreased from its value IT, the first peak in dash-dotted curve in Fig. [1] will move to 
lower frequency and so increases significantly the value of longitudinal conductivity at 
Q = 200cm~^. Its peak will cross this reference energy for 1/^B{T) ~ 1.47 and this 
produces the first peak in the long dashed (red) curve of Fig. [61 As 5 is decreased 
further below one Tesla, the second peak in the dash-dotted curve of Fig. [1] will also 
move through f2 = 200cm~^. This occurs for ~ 3.5, where a second lower 

intensity peak is seen in the long dashed (red) curve of Fig. [61 The other curves of Fig. [61 
can be traced out from similar considerations based on Fig. [H with attention paid to the 
evolution of these curves with changing value of B. 
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Figure 6. (Colour online) Real part of the longitudinal conductivity, Re<Jxx{^) at 
three frequencies, long dashed 200cm^^, dash-dotted 400cm^^, solid 800cm~^ as a 
function of the inverse of the square root of the magnetic field, with B in T. 

The other parameters are T = lOK, F = 15K and /i = 50K. 



3.2. Spectral weight 

An interesting quantity to consider is the optical spectral weight that falls between 
r2 = and Q = Qm with Q^n a variable upper limit in the integral 

W{n^)= / dnRea^^in). (26) 



For the case T = r = A = B = Owe have 



= —^\fi\6{Q) + f y - l/iM . (27) 

With Nf = 2 this leads to 

l^(a^) = -(7^|/i| + ^(^]„/2-|/i|)7^(^]„/2-|/i|))~-— ^ for fi^ > (28) 
h hi 

For finite 5, assuming for simplicity Mq < fi < Mi we get from Eq. fl23|) 

2c ^ Af„ + Af„+i ' 

Accordingly for A = 

e^NfV^p\eB\ ^ OiSl^ - M„ - M„+i) 



2c ^ M„ + M„+i 

n=0 



^ ^ + V ^ + 1 



Magneto- optical conductivity in Graphene 



18 



For large values the maximum that contributes to (150]) can be estimated from 



fi„ = Mat + Mn+i ~ 2MAr and so 

N 



8h\eB\v'j,/c' 



But 



and so 



v/iVTT ^ for > 1, 



h 2 



(31) 



(32) 



which agrees with Eq. (1251) as we would expect. We note that as we have stated, the 
area under each line goes as \fB, but the sum over the lines up to VLm gives another 
factor ~ 1/\/B which means that W{Vlrn) is independent of B. In Fig. [7] we show 
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Figure 7. (Colour online) Variation of the optical sum W{Q.m) (multiplied by h/e^) 
in cm^^ as a function of Vim in cm~^ for three different cases. Long dashed /i = 50K, 
S = OT and T = r = OK (see Eq. dash-dotted ^ = 50K, B = 0.2T, T = 0.5K 

and r = 4K; solid /i = 250K, B ^ 0.2T, T = 0.5K and T = 4K. 

numerical results for {h/e^)W{VLm) in units of cm~^ as a function of VLm is cm~^. The 
long dashed (red) curve is computed on the base of the first equality in Eq. (|28ll for 
/X = 50K, T = OK, and = OT and is given for comparison with the two other curves. 
These curves are obtained by numerical integration of Eq. ( 126|) with Yieaxxi^) given by 
Eq. ([021) for the case T = 0.5K, F = 4K, A = and finite B = 0.2T. The dash-dotted 
(black) curve is for chemical potential /i = 50K, while solid (blue) is for fi = 250K. 
For the first case fi is below Mi at 131.5cm~^ and for the second is above. At small 
values of both curves start at zero and remain small until Qm goes through the 
first Landau peak in Fig. [1] at which point it rises sharply and subsequently exhibits a 
plateau before showing a next sharp rise. For the dash-dotted curve (black) the first 
rise is at 131.5 cm~^ (interband line), but for the solid (blue) curve it is at the intraband 
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line (54.6 cm~^). As Qm increases to include several peaks in a^xi^) the curves start to 
follow more closely the B = OT result [long dashed (red) curve]. This behaviour can be 
understood better from a consideration of Eq. ( l30l) . The magnetic field mainly readjusts 
the available optical spectral weight among the Landau levels below ^2^- 

To be specific we have considered explicitly in this section only the case Mq < fi < 
Ml. For fi E]Miy, M]y+i[ we can show that the missing spectral weight in the lines 
n < is all to be found in the single intraband line at Mjy+i — M^r. From Eq. fl23l) 
and noting Eqs. (HM and (l25l) as T — > 0, the optical spectral weight lost in units of 
(eV/i)Mi(7r/2) is 

^'2 1 1 

The first term in the left hand side is the spectral weight from all the lines that have 
completely disappeared from n = to — 1. The second term is from the reduction in 
intensity by a factor of 1/2 of the line at n = A^. The quantity in the right hand side is 
the optical weight of the intraband line which has picked up all of the lost intensity. 



3.3. Hall conductivity 

Next we consider the absorptive part of the transverse Hall conductivity. The 
calculations are based on the equation (fTOj) . In Fig. [8] we show results for lmaxy{fl) in 
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Figure 8. (Colour online) The imaginary part of the Hall conductivity, lm<Txy{i^) in 
units of e'^/h as a function of frequency fl in cm"^. The three cases are for ^ = 50K 
(long dashed), /i — 510K (dash-dotted) and fi = 660K (solid). The other parameters 
are S = IT, T = lOK and T = 15K. 

units of e'^/h as a function of Q in cm^^. In all cases considered temperature T = lOK, 
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impurity scattering rate T = 15K and the excitonic gap A = 0. The magnetic field 
B = IT and three values of chemical potential are considered. The long dashed (red) 
curve is for /i = 50K below the energy of the n = 1 Landau level, dash-dotted (black) 
/i = 510K between n = 1 and n = 2, and solid (blue) fi = 660K between n = 2 and 
n = 3 as in Fig. [H We note a single peak in the long dashed (red) curve at Mi in 
contrast to two in the other two curves. For the dash-dotted (black) curve the peaks 
are at M2 — Mi and Mi + M2, respectively, and for the solid (blue) curve they are at 
M3 — M2 and M2 + M3. These features can be easily understood from Eq. ( ITOl) for 
a^yi^) when we take its imaginary part in the limit A = 0, F which is for i7 > 0, 
fi>0 



2 

Ima^yin) = -^M^^ J2[nF{Mn) - np{Mn+i) - np{-Mn+i) + ^^(-M^)] (34) 

n=0 

'5(^] - M„+i + M„) , 6{n-Mn+i~Mn 



h ' 2 

n=0 

M„+i + M„) 6{n - Mn+i ~ 
M„+i - M„ M„+i + M„ J ■ ^ ' 

Taking the limit of zero temperature, T = in Eq. only the first two thermal 
factors in the square bracket survive. For /i g]Mo, Mi[ we get 

Ima.,(fi) = -^MiV^^^^^, (36) 
while for /i E.]Mn, MAr+i[ with > we get 



, ,2 7r \5{n - Mn+i - Mn) ^ 6{VL - M^+i + M^j 



.(37) 



Mtv+i + Mat Mn+i - Mn 

Thus for /i e]Mo, Mi[ the transverse Hall conductivity exhibits a single peak at = Mi 
with weight 2/Mi, in units of {e^ /h){'K /2)Ml and for /i e]Mi^,MN+i[ with > 
it has two peaks at M^r+i -|- Mat and Mat+i — with weight l/(MAr+i + Mat) and 
l/(MAr+i — Mat), respectively. Thus the case n g]Mo,Mi[ is different from all others. 
This distinguishes Dirac from a classical Landau level quantization, for which, all cases 
would have two peaks. Also the area under the peaks in ImaxyiP) is ~ \fB in analogy 
to what we found for Yieaxxi^)- 

Next we consider the effect of a finite excitonic gap on the imaginary part of the Hall 
conductivity. In this case even forF— >0,f2>0,/i>0 equation (fTOj) is slightly more 
complicated than ( |34l) because of the additional factors 1± A^/M„M„+i which are only 
equal to 1 when A = 0. Nevertheless to understand physically the results given in Fig. [9] 
this extra complication is not needed and we can use Eq. flMl) as a guide. What is most 
important is to look at the thermal factors. We will consider only the case when the 
chemical potential falls between n = and n = 1 Landau level in energy. The relevant 
thermal factor is ufIMq) — nplMi). All four curves in Fig. [9] have n = 150K. The 
long dashed (red) curve is for A = OK and is included for reference. It shows a single 
peak (in —lmaxy{^)) at 294cm^^ as we know from our previous discussion. For the 
dash-dotted (black) curve A = lOOK which falls below the value of chemical potential 
(see Fig. [2] (b)). In this case the thermal factors in Eq. at T = give 1 and the 
delta functions correspond to = ^A^ + — A and = a/A^ + f2f + A. Thus, the 
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Figure 9. (Colour online) The imaginary part of the Hall conductivity, Im axy{^) in 
units of e'^/h vs frequency fl in cm"-'^ for temperature T = lOK, F = 15K, B = IT 
and chemical potential /i — 150K for 4 values of the excitonic gap A. Long dashed 
A = OK, dash-dotted A = lOOK, solid A = 250K, short dashed A = 200K. 

peak in the long dashed curve has spht into two. Further such peak has approximately, 
but not exactly the same optical spectral weight as that under the single peak of long 
dashed hne. This arises because of the weighting factors 1 ± A^/M„M„+i of Eq. (iTOll 
not shown explicitly in Eq. flM|) . For A = /i [solid (blue) curve] the two peaks remain 
at Q = a/A^ + i7f ± A (note the small shift which corresponds to the slightly different 
values of the gap between solid and short-long dashed curve). However, the thermal 
factor Up (A) - upiA + ^A^TTf) now equals 1/2 rather than 1 for the two previous 
cases, so that this feature on its own reduces the optical spectral weight of these peaks 
by half. Finally for the short dashed (green) curve the gap A is larger than is the 
chemical potential and the thermal factor nF{Mo) — upiMi) is zero, so that no peak is 
seen. 

In Fig. [To] we show results for the change in Ima^yi^) in units of e'^/h as a function 
of /i in Kelvin at fixed optical frequency. This figure is the analog of Fig. [5l The same 
parameters T = lOK, F = 15K, B = IT and A = are chosen as well as Q, namely, 
long dashed (red) curve, Q = Q2 — = 122cm^^. The absorptive Hall conductivity 
for this curve is near zero until the first Landau level energy 420K is crossed, where 
it drops below —17e^/h after which it shows a plateau till the next level is crossed at 
= 595K, where the second step is seen, etc. The dash-dotted (black) curve is for 
i7 = fii = 294cm~^. In this case the first plateau is at ~ —lAe^/h till the energy of 
the n = 1 Landau level is crossed in which case it drops to near zero value. The short 
dashed (green) curve is for VL = 0.95Qi = 280cm~^ and follows the dash-dotted (black) 
curve except the plateau is at ~ —lOe'^/h as we expect. Finally the solid (blue) curve 
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Figure 10. (Colour online) The imaginary part of the Hall conductivity, lmaxy{^) in 
units of e^/ft. as a function chemical potential /i in K. Four frequencies are considered, 
long dashed ^ fl2 — = 122cm^^, dash-dotted SI = fii = 294cm~^, solid 
n = ^2 + ^1 = TlOcm^i and short dashed Q = 0.95f^i = 280cm"^ The other 
parameters are B = IT, T = lOK and T = 15K. 



is for ^2 = ^21 + ^2 = 710cm^^. It starts at zero till /i crosses 420K at which point it 
shows a step down, remains nearly constant and finally steps back to near zero value at 
/i = 595K as expected from consideration of Fig. [HI 

The final Fig. [TT] shows results for the variation of lmaxy{^) in units of e^/h vs 
the inverse square root of the magnetic field, 1 / ^JB{T). Here temperature T = lOK, 
impurity scattering rate F = 15K and the chemical potential /i = 50K with excitonic 
gap A = 0. Three optical frequencies are chosen: long dashed (red) curve VL = 160cm~^, 
dash-dotted (black) curve Q = 240cm~^ and solid (blue) curve fl = 400cm~^. These 
curves are similar to those of Fig. [6] for the longitudinal case, but here exhibit only a 
single peak as is expected from the curves shown in Fig. [81 

4. Discussion 

In this paper we have extended the calculation of the ac conductivity of graphene 
presented in Ref. [H] in several ways. One emphasis has been on presenting results 
for different values of the chemical potential, /i. In a field-effect device, can be varied 
within the range of a few thousands of Kelvin by changing the gate voltage [H [2], [3l [T7]. 
This implies that for fields of the order of IT, /i can be made to sweep over several 
Landau levels. Recently [18] infrared spectroscopy has been successfully applied to study 
FET devices based on poly(3-hexylthiophene). We find a rich pattern of behaviour for 
the fl dependence of the real (absorptive) part of the diagonal conductivity, Recr^^^Q, T) 
as a function of in a fixed external magnetic field oriented perpendicular to the 
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Figure 11. (Colour online) The imaginary part of the Hall conductivity, Imaxyi^) 
in units of e^/h as a function of the inverse of the square root of the magnetic field, 
1/VB with B in T. Three frequencies are considered, long dashed 160cm^^, dash- 
dotted 240cm^^, solid 400cm~^. The other parameters are T = lOK, F = 15K and 
^ = 50K. 



graphene sheet. Peaks are seen in axx{^,T) corresponding to the possible transition 
between Landau levels from the hole to particle band (interband) or within a given cone 
(intraband). For /i between the n = and n = 1 levels there are lines at f2 = Mi, 
Ml + M2, M2 + M3 (interband), etc., where Mj is the position in energy of the i-th 
Landau level. There is no peak below Mi. The relative optical spectral weight between 
the various levels decay approximately as with fl evaluated at the center of each 
peak. When the chemical potential falls between the energy of the n = 1 and n = 2 
level the peak at f2 = Mi fades into the background and a new line appears at the 
lower frequency M2 — Mi (intraband). In addition, the intensity of the line at Mi + M2 
drops to half its value, while all other lines remain the same. As the chemical potential 
crosses higher and higher Landau level, say falls between M^r and M^r+i, the low energy 
peak has shifted to Mat+i — Mjv (intraband) after which all peaks previously seen have 
disappeared into background except for the one at Mjy+i+Miy which has half its previous 
intensity. Again all peaks above this energy remain unaltered. The peculiar behaviour 
of the peak at = Mi which is either present with full intensity or completely absent is 
the hallmark of Dirac as oppose to Schrodinger behaviour. This peak can be classified 
either as inter or intra when the excitonic gap A = 0, an ambiguity that is lifted for 
A 7^ 0. For Landau levels based on the Schrodinger equation the interband lines are 
evenly spaced while the intraband line is fixed in energy, and all interband lines first 
half their intensity before disappearing entirely as fi is increased. 

We have found that the transition from one configuration of absorption lines to 
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another occurs for a small change in chemical potential near a given Landau level 
energy with the scale for the incremental change in /i set by temperature and/or level 
broadening F. Away from these special values of /z the curves do not change significantly. 

In anticipation of experiments [32] we also provide scans of the the behaviour of 
^QO'xxi^jT) at fixed Q as a. function of chemical potential or of magnetic field. The 
pattern of behaviour found is traced back to that just described for Re a^xi^, T) vs fl at 
different values of field and chemical potential. Parallel results for the absorptive part 
of the transverse Hall conductivity are also presented. 

The possibility that an excitonic gap may open in graphene under high magnetic 
field has been considered by many authors [IHl EDI EH Elj and may even have been 
observed in recent experiments [23j. We have considered its effect on the absorption 
peak seen in both diagonal and Hall conductivity. Our specific predictions are that, 
a given peak can split into two, can disappear entirely or can simply shift to higher 
energy without splitting depending on the value of the chemical potential. In some 
circumstances the optical spectral weight under the splitted peak can differ from each 
other by a factor of order 3. 

Many of the results shown in the figures were obtained on the basis of general 
formulas for the conductivity, but it was found that simplified versions which can be 
more easily used to interpret experimental results are surprisingly accurate in the cases 
considered. They involve sums over Landau level indices of Lorentzian forms. Recently 
Li et al. [22] and Sadowski et al. [29] used related forms to analyze their results. 

Early magnetoreflectance data [33] in magnetic fields in the range 1 to 10 T carried 
out in graphite assigned lines to the if-point and some of these lines were found to follow 
a square root of B law as expected for graphene. Recent data of Ref. [22] at higher fields 
up to 20 T, however, found a conventional linear in B behaviour. Very recent infrared 
transmission data [2^ in ultrathin epitaxial graphite [301 [3T] in fields up to 4 T do 
show \fB law and these authors took this to be evidence that their carbon sheets are 
sufficiently decoupled to behave like graphene. However, this interpretation also requires 
that they assume that various sheets carry different charges, i.e. have different values of 
chemical potential. Nevertheless, the main lines seen in these data were assigned to the 
first three Dirac interband transitions and the (0, 1) intraband transition on the basis 
of their position in energy. They were also seen to vary in optical spectral weight as 
the square root of the magnetic field in good agreement with the findings in this work. 
Many other detailed predictions made in this paper have yet to be verified and should 
help in firmly establishing the special characteristic of quasiparticles in graphene. 
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Appendix A. Calculation of o'±{Q) in magneto-optical Lorentzian model 



We consider the complex conductivity (T±{Q) = o-^^^Q) ± ia^yi^). Substituting the 
spectral function ([3]) in Eq. ([2]) after evaluating tr and integrating over momentum k 
[see Appendix A of Ref. [8]] one obtains 

e^NfvlleB 



Aire fii 



" rf^rf^'!^£(^^^!L^£M [^,(u;,uj') T sgn(ei?)^2(^,^')] , (A.l) 



where the functions ipi{u!,uj') and ip2{^^,^^') are: 



n,m=0 



(A.2) 



with 
and 



A2 



+ 1 + 



An{uj) 



{An{uj)B^m{uj') + Bn{uj)A„,{uj')) 



Bnito) 



7r[(^-M„)2 + r2]' 
One can easily check that 

Using this symmetry one can rewrite Eq. (lA.ip in the form 



7r[(a; + M„)2 + r2] 



(A.3) 
(A.4) 
(A.5) 



e^Nfvl\eB\ 
Ancfli 

± sgn(ei?) 



1 



+ 



Lu ~ Lu' + Q + iO u — u' — Q — iO 
1 



uj ~ to' + + iO uj — to' — ri — iO 
Now using Eq. (]A.2I) we obtain 

e^Nfvl\eB 



(A.6) 



AncQi 



X 



n,m=0 J~OD 
1 1 



+ 



d= sgn(ei?) 



uj — Lu' + Q + iO u — uj' — Q — iO 
1 1 



cij — tu' + + zO Lu — u' — Q — iO 
For a typical integral over u' we have 



('^ra,m— 1 ~l~ ^m,n—l) 

('^n,m— 1 '^m.n— ly 



(A.7) 



duj'i:n,m{^^,Uj') 

LU -u' + il + iO 



AJuj) 



+ 



Bn{u:) 



+ 1 + 



— Mm + i7 + iVjn + Mm + f2 + iVm 

An{u) 



+ 



C(j — Mm + + + Mm + f2 + ir 



(A.8) 
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If the thermal factor n^;' is absent, further integration over cu would give an exact result: 



dujduj 



u; - cu' + Q + iO 



(A.9) 



1 - 



1 



+ 



1 



+ 1 + 



1 1 



MnMmJ \-Mn - M„, + Q + l{T ^ + Mn + Mrn + ^ + OV ^ + Tm)j' 

In case we integrate with a smooth function np iuj) we can approximately write 



/CX) pC 
duj Up {(jj) 
-oo J — c 



1 - 



A2 



00 -uj' + Q + iO ' 

Up (M„) 



(A.IO) 



+ 



+ 1 + 



M„M^ ; V M„ - + i7 + i{Vn + r^) -M„ + + 1] + i(r„ + r^) 



+ 



,-M„ - + Q + i(r„ + rj M„ + M„ + Q + i(r„ + rj 

Similarly, taking the complex conjugate and changing Q — > — Q, we have 

dljj'i)n,m{^,^') 



/oo /-c 
dujnp{uj) 
-oo J —c 



1 - 

1 



MnMm 



00 — u' — Q — iO 

np (Mn) 



(A.ll) 



np{-Mn) 



Mn-Mm-n- i(r„ + r^) -M„ + - n - i(r„ + r^) 

n^(-M„) , np{M„) 



-Mn-Mm-n- i(r„ + r^) m„ + - n - i(r„ + r„)y ' 



Hence we arrive at 

e^Nfvl\eB\ 



a±(Q) = 



X 



E 

n=0 



Attc Qi 
A^ 



(A.12) 



X 



[[rif (M„) - np (M„+i)] + [n^ (-M„+i) - (-M„)]] 

1 1 

- + 



M„ - M„+i - 1] - ^(r„ + r„+i: 



_Mn - M^+i + n + 2(r„ + 

J- ji+i 

A2 \ 

[[np (-M„+i) - rif (M„)] + [np (-M„) - (M„+i)]] 



- 1 + 



X 



+ 



±sgn(ei?) 



M„ + M„+i - Q - i(r„ + r„+i) M„ + M„+i + Q + i(r„ + r„+i) 

A2 \ 

[[np (M„) - up (Mn+i)] - [up i-Mn+i) - up i-Mn)]] 



X 



1 



1 



Mn - + n + i(r„ + r„+i) m„ - m^+i - n - i(r„ + r„+i) 



+ 1 + 



X 



) [[np (-M„+i) - np (M„)] - [n^ (-M„) - np (M„+i)]] 



1 



M„ + M„+i - (7 - i(r„ + r„+i) M„ + M„+i + n + i(r„ + r„+i; 
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where only the sum over n remained. 

We verified that Keaxxi^) computed for r„(ci;) = const from the more approximate 
Eq. ( ]A.12[) agrees quantitatively with the results obtained from a full Eq. ([71). This 
agreement is best for the resonance peaks and only small deviations are seen for 
r2 ~ 0. Nevertheless, Eq. (1A.12I) has a few drawbacks due to approximations made 
in Eqs. f lXlQ]) . flXTTj) . 

In particular, the Drude form cannot be recovered in i? ^ limit, while it can be 
obtained [HI [IE] from an exact representation ([7]). Moreover, the imaginary parts of the 
diagonal conductivity, Imaxxi^) ^^e Hall conductivity, lmaxy{^) are divergent in 
the limit ^ and do not satisfy Kramers-Kronig relations with the corresponding 
real parts found from Eq. flA.121) . To correct these problems we move the term 1/Q 
under the sum replacing it by its value at the pole of the corresponding denominator in 



Eq. (EH and arrive at Eqs. (El) and ([10]). 
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